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Abstract 



Qh! In this paper we investigate a possibility to perform shape-sensitive 



HBT analysis in transversal plane. We study presently used Side-Out 
decomposition of the correlation function and we do not find it to be 
very sensitive for the detection of a shape of source distribution in 
transversal plane. Therefore we suggest to use X-Y decomposition of 
the correlation function for the HBT analysis of non-central heavy ion 
. collisions. We test the method suggested on artificial and non-central 

I Pb-Pb 158GeV/n events produced by CASCUS event generator. Fur- 

■ ther possible enhancement of the method is sketched. 
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1. Introduction 

HBT interferometry analysis has become a powerful method in present day 
heavy ion collision (HIC) experiments. Nice theoretical and experimental work 



;h ' has been done in the field of HBT analysis of central HIC events. Long-Side- 

^ Out decomposition of correlation function enabled us to separate space and 

time properties of source distribution of bosons p{x,t). Even more sophisti- 
cated models of emission functions containing fiow and velocity parameters 
have been analyzed in various coordinate systems fl]]. 

On the other hand HBT analysis is still apparently not fully accomplished 
technique. For example Gamov factor correction which was expected to de- 
scribe infiuence of coulomb final state interactions on the correlation func- 
tion was found to be in disagreement with recent experimental results 0. 

Recently HBT analysis of non-central HIC events has been performed on 
data taken at SPS [Q] and AGS [Q. Extracted radius parameters were found 
to be dependent on centrality of collision and also on azimuthal orientation of 
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the momentum of pions in respect to the reaction plane [Q. For this analysis of 
non-central events a methods originally developed for cylindrically symmetric 
systems were used. 

We think that a more comprehensive information about the source distri- 
bution p{x, t) can be obtained. The purpose of this paper is to introduce an 
HBT analysis method which is suitable for the study of geometrical properties 
of source distributions e.g. a shape or structure of p{x). 

This paper is organized as follows: In Section 2 we study Side-Out decom- 
position of correlation function and we find it to be unsensitive to the shape of 
the source distribution in transversal plane. In Section 3 we introduce X-Y de- 
composition of correlation function and we show it to be suitable for the study 
of the geometrical properties of source distribution. A possibility to perform 
azimuthally sensitive HBT analysis of non-central events in the case when ori- 
entation of impact parameter b is not known is investigated in Section 4. In 
subsequent sections we study influence of x — p correlation of source distribu- 
tion p{x, t) on transversal HBT analysis and also possible further enhancement 
of the HBT interferometry method. 

2. Side-Out Decomposition of Correlation Function 

In this section we shall recall advantages of the Side-Out decomposition of 
correlation function. We start with simplified expression for the correlation 
function produced hj x — p non-correlated source p{x, t). Assuming that prob- 
ability of emission of a single boson does not depend on momentum and that 
emission of bosons is incoherent we write: 



C(Ap) = / / Sxi(fX2 p(xi,ti).p(f2,t2)|^i2r dtidt2 (1) 

where 



= 1 + cos(Af ■ Ap - At ■ A^) (2) 
By simple kinematic relations we can prove |Q: 

/\E = Ap- Pout ; Pout = , Ap = pi-p2 (3) 

Thus decomposing scalar product in Eq. (2) into Out and Side orthogonal com- 
ponents we can separate time information about the emission source from the 
Side component of correlation function: 

1 + COs{Apside ■ Ax side + Apout " [AXout " At ■ Pout]) (4) 

This Side-Out decomposition of correlation function allows to study time 
properties of source p{x,t) as a differnce of Rout^Rside extracted parameters. 
One of the main messages contained in this paper is: 
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Side- Out decomposition of correlation function destroys information 
about possible transversal shape of source distribution p{x, t) . 



Also for the sake of simplicity we concentrate on the study of transversal 
properties of source distribution p(x, t) in this paper. For this purpose we shall 
investigate properties of transversal correlation function constructed from pairs 
with very small longitudinal component of relative momentum \Api\ -C |Apt|. 
Then we can neglect longitudinal term in scalar product Ap ■ Ax in equation 
(2): 

Ap - Ax = Api ■ Axi + Apt ■ Axt ~ Apt ■ Axt (5) 
and for the transversal correlation function we write: 

C'^{Apt)= [ [ d'^xid'^X2\p\xi,t).p\x2,t)[l + cos{ApfAxt-At-AE)]\dtidt2 



(6) 

Source distribution p^{x) in (4) is a projection of three-dimensional source 
distribution p{x) onto the transversal plane: p^{xt) = J p{xt, Xi)dxi. 

On Fig.l we show transversal correlation function in Side-Out decomposi- 
tion for artificially asymmetric S-Pb events. Momenta of bosons were gener- 
ated in agreement with experimentally measured pt and rapidity distributions 
of pions in central S-Pb 200GeV/n collisions in the same way as it was done 
in computer simulation 0. Space distribution of emission points of pions in 
transversal plane was changed artificially as it is shown on Fig. la. 




Fig.l Side - Out Correlation function for asymmetric source in transver- 
sal plane. Pions are assumed to be emitted at the same time (At = 0). 

Transversal correlation function C^(Ap*j^g, Ajo^^J does not exhibit any az- 
imuthal asymmetry as a consequence of asymmetrical shape of source distri- 
bution in Fig. la - it is azimuthally non-sensitive. In the next section we show 
that another decomposition of correlation function is sensitive to the shape 
properties of p{x,t). Here we shall try to explain what is a mechanizm of 
destroying the shape information in Side-Out decomposition analysis. 
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a) b) 

Fig. 2 Two-point source emitting bosons. Relative momentum may 
be of different orientation than it is shown here. 

Let us study a very simple case - static two-point source emitting pairs of 
bosons simultaneously (Fig.2). If bosons are emitted as it is shown on Fig.2a 
then separation of emission points in Out and Side directions is: 

Axout = ; Axside = D (7) 
On Fig. 2b the orientation of sum pi + p2 leads to a totally different situation: 

AXout = D ; AXside = (8) 

Since the orientation of vector pout = 'eIX% random values of Axout, ^^side 
are randomly distributed but fulfilling condition Ax^^^+Ax^j^g = D^. It means 
that source behaves as a rotating in Side-Out transversal HBT analysis. 

In the next section we introduce X-Y approach to HBT which does not 
exhibit this feature and we find it to be sensitive to a non-trivial geometrical 
structure of source distribution in transversal plane. 

3. X-Y Decomposition of Transversal Correlation Function 

Let us imagine that orientation of impact parameter b is known with rather 
high precision for each non-central HIC event. In this case we can rotate 
measured momenta of particles in all events to have the same orientation of 
impact parameter. 

Because of asymmetrical overlapping region of colliding ions in non-central 
events the shape of source distribution is asymmetrical in transversal plane. 
This shape information would be destroyed if Side-Out type of HBT analysis 
without any additional cuts was applied^. 

On Fig. 3 we show transversal correlation function in X-Y decomposition 
for the same source distribution as shown on Fig. la. Correlation function 
C'^{Aprc, Apy) is clearly asymmetrical and this reflects the shape of original 
source distribution p{x,t). 

^See Section 5 for a possibility to apply directional cuts in momentum space. 
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b) C) 
Fig.3 X-Y transversal HBT analysis. Fig.Sc) is a contourplot of Fig.Sb). 
See also Fig.l for comparison. 

In the first approximation we suggest to use following gaussian parametriza- 
tion of the correlation function: 

C^(Ap,, Apy) = 1 + x.e-'^pl-Rl/^-^pl-Rl/^ (9) 
which corresponds to the static (sharp freeze-out) source distribution 

= c.e-^'/^'-^'/^? (10) 

Results of this simple X-Y decomposition HBT analysis for Pb-|-Pb non- 
central 158 GeV/n HIC events generated by CASCUS generator are presented 
in next section. Here we shall discuss the question of lifetime of source which 
was the main reason for Side- Out decomposition of C{Ap) as it is described 
in the preceeding section. 

Scalar product in the argument of cos in equation (2) in X-Y coordinate 
system (Ap^ 0) is: 



cos(Aa; • Ap^ + Ay ■ Apy - At ■ AE) 



It means that time information contained in p{x, t) is mixed in the same way 
with and Ry parameters of the source extracted by fit of the experimentally 
measured correlation function to equation (9). Consequently if source p{x, t) is 
asymmetrical in X-Y coordinate system then size parameters R^, Ry extracted 
by the fit should still exhibit asymmetry R^ ^ Ry. 

Moreover if we rely on the argument of cos in equation (2) then infiuence 
time dependet factor At ■ AE on our C^[Apx, Apy) formally vanishes if only 
pairs fulfilling the following condition 



AE — El — E2 







\Pi\ 



\P2\ 



(12) 



are selected for HBT analysis. Condition (12) together with our previous 
requirement Ap^ ~ still leaves enough freedom in the distribution of pairs in 
relative momentum space {Ap^, Apy). 
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We finish this section by schematic diagram depicting advantages of the 
above described Side-Out and X-Y HBT methods. 

Size R 
Lifetime t 
Flow T| 



Size R 
Shape • 
Structure © 

Fig. 4 Advantages of Side-Out and X-Y HBT analysis. 

Side-Out decomposition of C{Ap) allows to study mainly lifetime r and flow 
T] parameters of source (besides total size R). X-Y decomposition is sensitive 
to the shape and possibly other non-trivial geometrical structure of source 
distribution p{x). 

4. Pb+Pb 158 GeV/n Non-central Events CASCUS Simulation 

Azimuthally sensitive X-Y approach described in preceeding section was 
tested using a computer simulation of Pb-|-Pb 158GeV/n non-central events 
|6| = 7fm. Main structure of the simulation is shown on Fig. 5. 

GETIMI»ACT ► 

R 

CASCAmm ftESCATTERIKG , 

► ► ► 

GENERATOR PROGRAM ^ 

T 

mmmmmmmmmii Immmmmmmm J 

Fig. 5 CASCUS generator setup used for tiie simulation of non-central 
Pb+Pb 158 GeV/n HIC events. 

Initial positions and momenta of pions created in non-central Pb-|-Pb col- 
lisions were produced by Cascading Generator Pions were evolved in time 
using enhanced version of rescattering program [^. Resulting distribution 
of last interactions of pions was used for HBT analysis. 

Events with random orientation of impact parameter \b\ = 7fm were ro- 
tated to have the same orientation of impact parameter b. Cascading genera- 
tor program was run at Czech Academy of Sciences (CAS) and rescattering 
program |^ was written at Comenius University in Slovakia (CUS). 

Correlation function was constructed in the way similar to that used in 
simulation but in X-Y coordinate system. Result was fitted to analytical 
expression (9). Parameters Rx,Ry obtained by the fit are summarized in the 
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following table. 



X-Y HBT 


p{x) 

UHiVoriiN AIj 


BEFORE 
ilii/bL'Ai lii/rtiiNVjr 


AFTER 
rtil/bUAl lii/rtiiNLj 


TOY 
EVENTS 


Sx = 16/m 
Sy = Afm 
Sx - Sy = 12/m 


Rx = 8.76 ± 0.02/m 
Ry = 2.35 ± 0.02/m 
AR = 8.4/m 


Rx = 8.12 ± 0.02/m 
Ry = 3.67 ± 0.02/m 
AR = 7.2 fm 


Pb + Pb 
EVENTS 


Sx = 11/m 
Sy = 6/m 
Sx - Sy = 9/m 


Rx = 5.28 ± 0.05/m 
Ry = 3.44 ± 0.05/m 
AR = 4.0/m 


Rx = 5.48 ± 0.05/m 
Ry = 4.64 ± 0.05/m 
AR = 2.9 fm 



Tab.l Results of X-Y HBT analysis. Rx,Ry parameters are extracted 
by fit to Eq.(9). Sx, Sy are size parameters of source distribution p{x). 

As a consequence of asymmetrical shape of source distribution p{x) ex- 
tracted parameters Rx,Ry are not of the same value what is demonstrated 
also by non-zero value of AR = ^R^ — Ry in the above table. 

Comparing results obtained by HBT analysis of the events before and after 
rescattering proces we conclude that rescattering process smears the shape of 
the original source distribution. On the other hand x — p correlation in the 
resulting source distribution p{x,p) after the rescattering process allows to 
perform another type of shape sensitive analysis based on directional cuts in 
momentum space (see next section). 

5. Influence of x — p Correlation 

Let us imagine that there is no if— p correlation in the distribution of pions 
before the rescattering process. It means we assume that pions are created in 
individual nucleon - nucleon collisions and that there is no collective behaviour 
of nucloens before the emission of pions. Then geometrical properties of the 
emission points distribution of pions p{x,p) = p{x).p{p) do not depend on cuts 
in momentum space. On the other hand rescattering process creates x — p 
correlation in last interactions distribution of pions p^\x,p) (see Fig. 5). 
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Fig.5 Emission points of pions (CASCUS Pb+Pb 160 GeV/n ; b = 7fm) 
in transversal plane. Distribution of Last Interactions depends on pt cut. 
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Distribution of last interactions depends on absolute value of transversal 
momentum. This explains dependence of extracted HBT radii parameters 
on pt which was measured experimentally 0]. Moreover rescattering process 
creates also a directional x—Hp type of correlation ||10|. This type of correlation 
can also be used for detection of the transversal shape of p^^(x, p) distribution. 
On Fig.6 we show p^^{x,p) distribution of pions for selected directional cuts 
in transversal momentum space. 
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Fig.6 Directional cuts in momentum space select pions from different 
regions of Last Interactions distribution. Distribution of pairs in relative 
momentum space (Apt) still allows two-dimensional HBT analysis^. 

Pions with different directions of transversal momentum come from differ- 
ent regions of If distribution p^\x) is asymmetrical in transversal 
plane (see Fig.6) then directional cuts in momentum space select spatial re- 
gions of different size what leads to different radii parameters extracted by 
HBT method. Azimuthally sensitive HBT analysis based on directional mo- 
mentum cuts was already performed on AGS data [ffl. 



•^Directional cut in transversal momentum space {px,Py) does not restrict distribution of 
pion pairs in relative momentum space (Ap^;, Apj,) significantly. 
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Common feature of the method based on x — np correlation and X-Y HBT 
is that orientation of impact parameter is determined for each event (e.g by 
transversal flow analysis). In the next section we sketch a possibility to per- 
form azimuthally sensitive HBT without information about impact parameter 
orientation. 



5. Can We See the Reaction Plane via HBT ? 

In experiment the orientation of impact parameter b is random for non- 
central collisions. Thus transversal shape of source distribution will not be 
visible by HBT if the correlation function is constructed as a sum of randomly 
oriented single event contributions. This situation is similar to transversal 
flow analysis |]ll| where event by event (EbyE) technique allows to determine 



a presence of azimuthal asymmetry |l! 

Here we sketch EbyE HBT technique similar to that utilized in transversal 
flow analysis [jl3|. Let us rewrite X-Y decomposition of transversal correlation 
function 

C^(Ap,, Apy) = 1 + e-^P^«^/2-'^P?4/2 (13) 
in spherical coordinate system: 

Ap^. = Apr ■ cos((/)) ; Apy = Apr ■ sin(0) (14) 
After substitution (14) into (13) we obtain: 

C^iApr, 0) = 1 + e-¥k+Ai?^-(2<^)] (15) 

where = Rl + ; AR^ = Rl — Ry. Thus non-zero value of AR param- 
eter is a signature of transversal asymmetry of the correlation function and 
consequently also of the source distribution. Random orientation of impact 
parameter means a simple random shift in azimuthal angle <p — > + A* in 
spherical coordinate system. Equation (15) can then be rewritten in the form: 

C^(Ap,, 0) = 1 + [^O+E^ cos(2^)+i^, sin(20)] (^g) 

where parameters = AR^ cos(A*) ; Ey = AR^ sin(A*) depend on A* 
which is different for each event. This expression for single-event correlation 
function is similar to the second order fourier decomposition of azimuthal dis- 
tribution R{(fi) used in transversal flow analysis |I3| : 



i?(0) =Xo + X2 cos(20) + y2 sin(20) (17) 



Non-zero positioned maximum in f 2 = yX2 + distribution constructed 
from many single event fit results is a signature of non-zero value of average V2 
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parameter [O . Based on this analogy we imply that non-zero positioned max- 



imum of /S.E? = \J El + E"^ distribution obtained by many single-event HBT 
fit procedures could be a signature of asymmetrical shape of p{x) distribution 
in transversal plane. 

In this EbyE approach the orientation of impact parameter is not deter- 
mined by other (e.g. flow) method. Therefore we dare to say that experi- 
mentally observed non-zero average value of Ai? parameter extracted by this 
method would mean observation of the reaction plane via HBT. 

We are aware that here we give a very rough description of this method. We 
have not touched other important topics e.g. the construction of single-event 
correlation function or final state interactions correction IT^ . 



There exists also a possibility to perform HBT search for the reaction plane 
based on x — rip correlation described in the preceeding section. This approach 
would be based on the study of (anti) correlation of the radii R^^, R^^ extracted 
for e.g. orthogonal directional cuts in momentum space. 

These topics are however beyond the scope of this paper now. In the next 
section we describe shortly a previously reported [|l^, [l^ possibility of further 
enhancement of HBT analysis. 



6. Inverse Transformation 

General principle of HBT is the following: We measure a result of interfer- 
ence phenomenon (correlation function C{Ap)) and by analysis of C{Ap) we 
reconstruct some information about a source emitting bosons. 

This principle is similar to holography methods where however full spatial 
information about source can be obtained. Full reconstruction is allowed by 
the interference of coherent radiation. In most of literature about HBT a state- 
ment for degenerated (constant) correlation function in the case of coherent 
emission of bosons can be found. However it is not excluded that for coherently 
emitted bosons (produced by e.g. Chiral Condensates |T^) a slightly different 
- holographic^ {^^) approach to HBT exists. Nevetherless experimental sit- 
uation is the main Veto device which decides the fate of theoretically existing 
methods. Therefore we shall deal only with uncoherently radiated bosons in 
this section. 



As it was shown in work [|15| for a static (sharp freeze-out) x—p uncorrelated 
source we can write: 

C{Ap,,Apy)= 1 1 SiD,,Dy)[l + cosiAp.,- AD, + ApyADy)]dD,dDy (18) 

where S{D,, Dy) is a probability distribution for emission of bosons at a rel- 
ative distance D = {D,,Dy). Inverse transformation from C{Ap,, Apy) to 

^For the inverse transformation based on Helmholtz-Kirchhoff theorem see IE0| . 
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S{Dx, Dy) in equation (18) exists. It is quite tempting to perform this proce- 
dure on our simulated CASCUS events. 

On Fig. 7 we show S{Dx, Dy) - symmetrized result of inverse transformation 
of correlation function C{Apt) for artificial S-Pb events (see Fig. la) in X-Y 
a) and Side-Out b) decomposition. Images c), d) are result of the same 
procedure performed for two-point source (see Fig.2). 





a) 



b) 




o 



c) 



d) 



Fig. 7 Results of inverse transformation of equation (18) for correla- 
tion function of toy S-Pb events in a) X-Y (Fig.Sb) and b) Side-Out 
(Fig. lb) decomposition. Images c), d) are obtained by the same pro- 
cedure for two-point source (Fig.2). Distance between two dots in X-Y 
decomposition Image c) is 6fm. All images have the same scale. 

Distance distribution S'(D)0 in X-Y decomposition reflects clearly shape of 
the original source distribution shown on Fig. 3a. This feature allows to study 
a transversal shape of the source distribution in non-central collisions. Image 
Fig. 7c) confirms that X-Y HBT method is sensitive also to a non-trivial struc- 
ture of p{x). This feature might allow to perform a search for non-statistical 
inhomogenities possibly generated by a phase transition in central heavy ion 
collisions 



19 



^See [fM ITbI for the explanation of symmetrized S{D) distribution. 
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Conclusions 



We have shown that X-Y decomposition HBT analysis is suitable for the 
investigation of geometrical properties of the source distribution. This opens a 
possibility to study shape of the volume radiating bosons in non-central heavy 
ion collisions and possibly also non-trivial structure (inhomogenities) in central 
heavy ion collisions. 
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